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We study effects of nonmagnetic impurities in a spin-1/2 frustrated triangular antiferromagnet 
with the aim of understanding the observed broadening of ^^C NMR lines in the organic spin liquid 
material k-(ET)2Cu2(CN)3. For high temperatures down to J/3, we calculate local susceptibility 
near a nonmagnetic impurity and near a grain boundary for the nearest neighbor Heisenberg model 
in high temperature series expansion. We find that the local susceptibility decays to the uniform 
one in few lattice spacings, and for a low density of impurities we would not be able to explain the 
line broadening present in the experiments already at elevated temperatures. At low temperatures, 
we assume a gapless spin liquid with a Fermi surface of spinous. We calculate the local susceptibility 
in the mean field and also go beyond the mean field by Gutzwiller projection. The zero temperature 
local susceptibility decays as a power law and oscillates at 2kF- As in the high temperature analysis 
we find that a low density of impurities is not able to explain the observed broadening of the lines. 
We are thus led to conclude that there is more disorder in the system. We find that a large density 
of point-like disorder gives broadening that is consistent with the experiment down to about 5K, 
but that below this temperature additional mechanism is likely needed. 

PACS numbers: 



I. INTRODUCTION 

Spin liquid phases are some of the most interesting 
phases known to exist theoretically. However they are 
hard to achieve experimentally because interactions usu- 
ally favor ordered phases. To achieve spin liquid we have 
to frustrate these interactions. Triangular lattice pro- 
vides a natural way to do this. For the nearest neighbor 
antiferromagnetic Heisenberg model the frustration is not 
strong enough and the ground state is ordered. However, 
the order is weak, and it is likely that in the presence 
of appropriate additional interactions, spin liquid phases 
arise. 

This work is motivated by the layered organic 
compound K-(ET)2Cu2(CN)3.ii^i^iii^i^ It contains ET 
molecules that pair up, each pair lies on sites of tri- 
angular lattice and has one electron less then the full 
filling. The material at ambient pressure is an insula- 
tor. Thus it is effectively a spin 1/2 antiferromagnet on 
the triangular lattice. While the Heisenberg exchange 
J w 250 K, the system shows no signs of ordering down 
to 32 mK making it a good candidate for the spin liquid. 
There are likely additional interactions among spins, es- 
pecially ring exchanges that are thought to be respon- 
sible for driving the system into the spin liquid i^i^iiSiii^ 
What makes the appearance of such interactions natural 
is that under moderate pressure the k-(ET)2Cu2(CN)3 
undergoes a transition to a superconductor at low tem- 
perature (and a metal at higher temperature), so there 
are significant virtual charge fluctuations present already 
in the insulator at ambient pressureJ^ Crudely, we can 
think of the system as a half-filled Hubbard model close 
to the Mott insulator - metal transition, and we can esti- 
mate that the ring exchange interactions in the effective 
spin model are strong, enough to destroy the magnetic 
orderii^iii An alternative explanation of the insulator 



in terms of inhomogeneous electron localization has also 
been suggestedj^i^ 

The spin liquid phase remains enigmatic. Thermody- 
namic measurements show many gapless excitations in 
this charge insulator - at least as many as in a metal. One 
appealing proposal that captures some of the observed 
phenomenology is a state with spinon Fermi surface i^i^iii^ 
Other scenarios have also been suggested fi^'i^'i^ii^ par- 
ticularly with the view towards low temperatures. 

We are specifically interested in the ^'^C NMR mea- 
surements of the Knight shifts^i^i^ and what we can learn 
from these for the material and the spin liquid. The mea- 
surements effectively give a histogram of local magnetic 
susceptibility and are therefore a good probe of the mag- 
netic properties. The experiment shows strong broad- 
ening of such histogram as one lowers the temperature. 
The width of the peak broadens by about a factor of 
40 as the temperature is lowered from 250 K to roughly 
1 K and saturates as the temperature is lowered further. 
The distribution of local susceptibilities can be produced 
by disorder, since the susceptibility can have various val- 
ues as a function of distance from say an impurity. It is 
hard to imagine other mechanism producing a distribu- 
tion (except spin glass, but no such behavior is observed). 
Therefore in this paper we investigate the effects of dis- 
order on the spin system. 

Unfortunately, not much is known about the impurities 
and their role in the insulating k-(ET)2Cu2(CN)3 at am- 
bient pressure. At high pressure of 0.8 GPa, the material 
is a relatively clean metal with kpl^ 50— 100 and observ- 
able Shubnikov-deHaas oscillations. It is believed that 
the Cu^+ impurity concentration is very low^ less than 
0.01%. There are quite possibly additional sources of dis- 
order such as different local environments coming from 
different conformations of the ET molecules , or 
from disorder in the insulating anion layers such as the 
disordered (CN)~ group i^i^iSii^ Analyzing the NMR ex- 
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periments can then provide some understanding of the 
disorder, its strength, and role in the insulator phase. 
Taking up the Mott insulator picture as one viable candi- 
date, where the insulator is primarily driven by electron- 
electron interactions, we set out to study models of non- 
magnetic disorder in a spin- 1/2 system on the triangular 
lattice. We study progressively different kinds of disorder 
and analyze what each predicts about the local suscepti- 
bilities in turn. 

This paper consists of two separate approaches: one, 
the high temperature series expansion, and the other, low 
temperature analysis assuming the system forms a spin 
liquid. High temperature series expansion is rather re- 
stricted by the range of temperatures and the types of 
models it can study, but for those models it gives exact 
results. We are able to calculate local susceptibilities of 
the nearest neighbor Heisenberg model on an arbitrary 
graph where all exchange couplings are the same. In 
particular, wc can study triangular lattice with a miss- 
ing site, with a boundary, or in the presence of a finite 
density of missing sites. The missing sites are models of 
nonmagnetic impurities, while the boundary is a model 
of grain boundary. We can go down to temperature of 
roughly J/3. In this range, the experiments already see 
broadening of the Knight shift distribution by about a 
factor of two and thus we can compare the calculated re- 
sults to the experiments. Of course, the real material has 
interactions beyond the nearest neighbor since it is a spin 
liquid, but rough estimate can be made, especially given 
that multi-spin exchanges are less important at high tem- 
peratures. 

At low temperatures we assume phenomenologically 
the system forms a spin liquid with Fermi sea of spinous 
(stabilized by additional interactions). We first analyze 
this in mean field where it reduces to free fcrmions hop- 
ping on the triangular lattice. The spin liquid can nat- 
urally accommodate nonmagnetic disorder in the spin 
model by the corresponding changes in the spinon hop- 
ping amplitudes. The full theory also contains a dynam- 
ical U(l) gauge field^i^i2ii25j26^ ^^lis is hard to an- 
alyze directly. Instead, to go beyond mean field we study 
wavefunctions obtained by Gutzwiller projection of the 
mean field states. In one dimension, this can capture 
the full theory, while in two dimensions this is only an 
approximation but a reasonable one and dealing directly 
with physical spins. 

Overall, we find that the local susceptibility decays 
rather quickly near an impurity and at small impurity 
densities such as 0.01% of Cu our results are very far from 
explaining the experimental data - they would produce 
very sharp histograms as most of the sites are in the bulk. 
In the spin liquid phase, the local susceptibility has an 
oscillatory 2kF component that decays with a power law 
envelope away from defects, so the impurities can be felt 
at larger distances, but the overall amplitude that we find 
is still small. 

We then studied the system at high temperature near 
a boundary and in the presence of a larger density of 



missing sites. We also studied the system at low temper- 
ature in the mean field near a boundary, in the presence 
of a larger density of missing sites, and in the presence 
of random disorder on bonds which is either uniform 
or localized at a fraction of bonds. From this analysis 
it appears that the most likely scenario is the case of 
relatively large density of point-like disorder, where the 
lincwidths broaden with lowering the temperature until 
the correlation length becomes comparable to the typi- 
cal distance between defects. A puzzling feature is that 
with such fixed disorder we cannot reproduce the ob- 
served strong temperature dependence of the NMR lines 
as the temperatures are lowered further. On the other 
hand, such models in the metallic phase where we consid- 
ered electrons with random on-site potentials match rea- 
sonably with the experiments under pressure, where the 
NMR linewidths remain unchanged with temperature.— 
It could be also that the effective strength of disorder 
increases as the temperature is lowered in the insulator, 
perhaps because of the vanishing screening of charged 
impurities. Better understanding of the disorder in the 
K-(ET)2Cu2(CN)3 system is clearly needed. 

A new triangular lattice spin liquid material 
EtMe3Sb[Pd(dmit)2]2, Ref. [2^, has rather similar phe- 
nomenology to that of the k-(ET)2Cu2(CN)3 and also 
appears to have significant NMR line broadening, which 
may thus be a common feature of gaplcss spin liquids. It 
would be interesting to compare both systems more. 

Finally we would like to mention that we have per- 
formed similar analysis on Kagome antiferromagnet ad- 
dressing the NMR line broadening in the candidate spin 
liquid material ZnCu3(OH)6Cl2.^'20i3i There, the disor- 
der is relatively well understood experimentally and is 
estimated to be about 5% vacancies. Our calculations^^ 
in this case compare sensibly with the experiment. 



II. SUMMARY OF THE EXPERIMENTAL 
INHOMOGENEOUS LINE BROADENING 

In what follows, wc calculate local susceptibility in spin 
models with nonmagnetic disorder. To set the stage, 
we summarize the main experimental findings in the 
form convenient for judging theoretical results. A di- 
rect comparison with the experiments is to look at the 
width of the local susceptibility histograms relative to 
the average susceptibility. In the model calculations the 
bulk susceptibility is roughly the location of the his- 
togram peak, while in the experimental plots we should 
also be aware of the chemical shifts. The referencing 
to the average susceptibility is justified since in the n- 
(ET)2Cu2(CN)3 this remains roughly unchanged around 
X = 5 • 10~^ emu/mol in a wide temperature range be- 
tween 300 K and 30 K and then decreases somewhat to 
a value around x = 3 • lO""* emu/mol at 1 K. Also, the 
bulk values can be quantitatively reproduced by suitable 
choices of the model parameters such as J in the high- 
temperature series studjsiiii or the spinon hopping am- 
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plitude in the spin liquid model (see Sec. lIVp . 

Refs. Ulili show the NMR lines plotted versus shifts 
from tetramethylsilane (TMS), in parts per million 
(ppm). For the more strongly coupled ^^C whose hy- 
perfine coupling constant is 0.21 Tesla/(//B dimer), the 
susceptibility of x = 5 • 10~^ emu/mol corresponds to 
SB/B = 1.8 • 10-4 = 180 ppm Knight shift. Thus 
a 20 ppm shift corresponds to roughly a 10% relative 
change in the susceptibility at higher temperatures and 
a somewhat larger relative change at lower temperatures. 
Reading from the experiments, the full width at half 
maximum (FWHM) of the line measured in such relative 
terms increases from about 10% at 300 K to about 20- 
30% at 50 K to about 50-60% at 10 K, and then increases 
steeply to about 300% at 1 K and saturates around this 
value at still lower temperatures. One can appreciate 
the dramatic broadening directly from the line shapes 
m Refs. illE where the distance of this line from the 
origin sets a natural scale since the Knight shift and the 
chemical shift are comparable. (Note that to see the inho- 
mogeneous Knight shifts over the intrinsic linewidths the 
experiments arc done in large fields of order 8 Tesla; this 
may be inducing more significant ground state modifica- 
tions, while in this work we focus on the linear response 
susceptibilities and their distributions.) 

We note that the broadening first sets in gradually 
coming from high temperatures. The line is already 20- 
30% broad at T = 50 K ^ J/5 , which we can hope to un- 
derstand quantitatively using reliable high-temperature 
series approach. It further broadens by about a factor of 
two to three in the region 50 K to 5-10 K where we ex- 
pect the spin liquid approach to become applicable. We 
now present these two studies. 



III. HIGH TEMPERATURE SERIES 
EXPANSION 

We consider spin 1/2 nearest neighbor anti- 
ferromagnetic Heisenbcrg model on triangular lattice. 
Local spin susceptibility at site i is given by 



Xloc(«) 



(1) 



where S^^^ = "Y^- 5J. We calculate x's in the high tem- 
perature series expansions in the presence of nonmag- 
netic impurities treated as missing sites (vacancies), and 
also near an open boundary which is a model for a grain 
boundary. 

The expansion is performed to the 12-th order in J (or 
13-th order in 1/T) using the linked cluster expansionist 
The outline of the procedure is as follows. We generate all 
abstract graphs up to desired size. Then we generate all 
subgraphs of these graphs, keeping track of the location 
of each subgraph in the graph. 

We calculate the local susceptibility for each graph at 
each point of the graph. Then we subtract all the sub- 
graphs of each graph as needed in the linked cluster ex- 



pansion to get the contribution this graph would give 
when embedded into lattice. The local susceptibility on 
any lattice can be calculated by creating all possible em- 
beddings of all graphs and adding their contributions at 
every site. In this general formulation the lattice does no 
need to be regular, it can be any connected graph. In 
particular this procedure applies to the triangular lattice 
with a missing site, with a finite density of missing sites, 
or with a boundary. In practice, for the single impurity 
case, we consider all the graphs containing the impurity 
and subtract their contribution from the uniform suscep- 
tibility. Similar procedure is used in the case with the 
boundary. In the end, we obtain exact l/T series expan- 
sion for the system with such disorder 

After obtaining the scries, we extend it beyond the 
radius of convergence using the method of Pade approxi- 
mants. We use [5,6], [5,7], [6,6], [6,7] and expand in vari- 
able 1/{T + a) where a is usually 0.08 as in Ref. [s^. De- 
pending on a one might get a pole in the expression, and 
hence divergence in susceptibility even at relatively large 
temperature. This usually happens say in one of the ap- 
proximants while the others still overlap. At low enough 
temperatures they start diverging and we take that as a 
point where the approximation stops being valid. Differ- 
ent values of a are tried, and sometimes it is possible to 
tune to a value where all the curves overlap completely to 
a much lower temperature, but that is a pathology, prob- 
ably indicating that the polynomials are all the same. For 
other values of a the curves usually start to diverge from 
each other at around the same temperature. 



A. Point Impurity and Nonzero Density of 
Impurities 

The coefficients of the susceptibility of the first five 
nearest neighbors near the impurity are in Table [B The 
corresponding local susceptibilities along with the uni- 
form susceptibility are plotted in Figure [TJ We sec that 
the local susceptibility decays to the uniform one in few 
lattice spacings. This is consistent with calculated very 
short correlation length in Ref. [H. The deviation of the 
near-neighbor local susceptibility from the uniform value 
reaches roughly 15% at T « J/3. 

We would like to see if the observed NMR lines can be 
explained from a finite density of such impurities. The 
prediction is obtained by plotting the histogram of sus- 
ceptibilities. First we note that at T ^ lOOK ~ J/3 
the experimental lines are spread by about ±10% which 
is roughly equal to the calculated deviation of the local 
susceptibility of the nearest neighbors from the uniform 
susceptibility. The experimental curves have a significant 
weight spread over this width, and so if the calculated 
curves are to explain them, the first few nearest neigh- 
bors of impurities would have to form a sizable fraction 
of the total number of sites. Ref. [l| suggests that the 
system contains about one Cu impurity in ten thousand 
which gives about 0.4% fraction for the nearest neighbors 
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TABLE I: Series coefficients an of x ~ ^ X/n=o A^+^{n+i)\ ( t)" ^'^^ susceptibility Xuniform of tlie pure triangular lattice and 
for the five closest neighbors Xi of nonmagnetic impurity as indicated in Figure [T] The coefficients for the uniform susceptibility 
agree with Ref. [H. 
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Xuniform 


Xi 


X2 


Xa 


X4 


X5 
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1 
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-12 


-10 


-12 


-12 


-12 


-12 


2 


144 


108 


132 


138 


144 


144 


3 


-1632 


-1248 


-1312 


-1400 


-1560 


-1608 


4 


18000 


15840 


13840 


13320 


14880 


16160 


5 


-254016 


-237024 


nor— i-7i-7r> 

-235776 


-213984 


-189168 


-192000 


6 


5472096 


4144000 


5539968 


5817504 


5084464 


4564560 


7 


-109168128 


-73210624 


-93128960 


-109647744 


-123994240 


-118354560 


8 


818042112 


1133266176 


222006528 


112173696 


913327488 


1312247808 


9 


17982044160 


-18170275840 


11644656640 


30128806400 


38868680960 


28664414720 


10 


778741928448 


581215033344 


1535178191360 


1512448745984 


328581324544 


-104688021504 


11 


-90462554542080 


-21239974981632 


-84715204509696 


-115649955864576 


-118987461639168 


-96786926315520 


12 


829570427172864 


215676565092352 


-788032311226368 


-332026092103680 


2149211723363328 


2738259718125568 
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FIG. 1: [Color online] The uniform susceptibility and the local 
susceptibilities at the first five inequivalent neighbors of non- 
magnetic impurity (missing site). The four curves for each Xi 
are the Fade approximants [5, 6], [5, 7], [6, 6], [6, 7]. x is in units 
of {giiBf/J- 



up to X5 in Fig. [T] and hence it is very far from explain- 
ing the experimental lines - it would predict very sharp 
histograms. 

One possible explanation of this discrepancy is the fact 
that the Heisenberg model is not entirely adequate be- 
cause it would eventually predict ordering at low tem- 
peratures, which is not observed in experiments, and so 
there are additional interactions. Indeed, Ref. H proposed 
that ring exchanges are important to stabilize the spin 
liquid phase. However, since these interactions are still 
short-range, the volume fraction of sites that are affected 
by vacancies at these temperatures is still small, so low 
impurity density cannot fit the observed data. 

A different more likely explanation is that there are 
more impurities or more disorder in the system. One 



possible source of disorder is from extended defects such 
as grain boundaries, and in the next section we con- 
sider the susceptibility near a boundary. Other possi- 
ble source is from the ethylene group disorder in the 
ET molecules which is thought to be important in a re- 
lated K-(ET)2Cu[N(CN)2]Br material -i^iiiii^iiS This can 
give rise to a large density of point perturbations which 
are mild but present everywhere. Other possible source 
is disorder in the (CN)~ groups in the insulating anion 
layers We currently do not know much about the 
presence and magnitude of such perturbations in the spin 
liquid K-(ET)2Cu2(CN)3 material. 

To simulate of a large density of point-like dis- 

order, we study the system in the presence of a 5% of 
missing sites. Realistic point-like disorder is probably 
of a different nature, but vacancies is all we can do in 
the systematic high temperature expansion. However we 
hope the basic features of the histogram would be sim- 
ilar. The result is shown in the Fig. [21 Due to a large 
number of diagram imbcddings we were able only to go 
to the 11th order in J (rather than 12-th as above) on a 
20 X 20 lattice and rehably only down to T 0.5 J. The 
error on the histogram values is roughly ±20%. Crudely, 
we sec two sets of peaks in Fig. [2J The one on the right 
is associated with the nearest neighbors of the impuri- 
ties while the one on the left with the rest of the sites. 
The finer features are associated with sites at different 
positions with respect to several impurities. 

As we have already mentioned, the nearest neighbor 
susceptibility is different from the the uniform one by 
about 10%, which is roughly similar to the broadening 
of the NMR line at T - 100 - 50 K. However the tem- 
perature dependence of xi ~ Xuniform is also weak, see 
Fig. [TJ In the case with very low density of impurities 
such as 0.01% Cu, besides having this value only at very 
few sites, it would not give broadening of the lines. In 
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the case of 5% of impurities, Fig. [2l this corresponds to 
the distance between the two sets of peaks not changing 
significantly with temperature. On the other hand the 
(set of) peaks themselves visibly broaden more. Remem- 
bering that here we plot histograms corresponding to 
the "ideal" system (i.e., with only Heisenberg exchanges 
and vacancies, and not including other interactions and 
sources of line widths), few percent of impurities can in- 
deed produce reasonably broadened lineshapes at these 
elevated temperatures. 
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FIG. 2: The histogram of local susceptibilities for a 20 x 20 
triangular lattice sample with 5% randomly placed vacancies, 
obtained from the high temperature series expansion to 11-th 
order in J. x is in units of [g^B)'^ /J- 



B. Line Impurity 

In this subsection we consider pure triangular lattice 
near a boundary and calculate susceptibility at various 
distances from the boundary. The coefficients of the uni- 
form susceptibility and of the local susceptibility at the 
first five closest inequivalent sites are in Table |TT1 The 
corresponding plots are in Figure [31 

We see that the ratio of the local susceptibility to the 
uniform one is somewhat larger for the first neighbor here 
than in the single vacancy case and more importantly 
the temperature dependence is stronger. Furthermore, 
the line impurity is an extended object, so a much larger 
number of sites is affected. We do not know how com- 
mon such grain boundaries are in the k-(ET)2Cu2(CN)3 
material to make predictions for the experiment. How- 
ever, as in the case with vacancies, we see that down to 
T ~ J/4 only the first few Xi near the boundary deviate 
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FIG. 3: [Color online] The uniform susceptibility and the lo- 
cal susceptibilities at the first five inequivalent neighbors of 
the boundary. The four curves for each x the Pade ap- 
proximants [5, 6], [5, 7], [6, 6], [6, 7]. x is in units of {gfis)^ /J. 

significantly from the bulk value. Thus to explain the 
experimental lines with such defects we would require a 
large density of them. 

IV. SPIN LIQUID WITH FERMIONIC SPINONS 

The k-(ET)2Cu2(CN)3 system does not order down 
to temperatures as low as 32 mK, but has a J « 
250 K. Thus it is a good candidate for spin liquid. 
Among SU(2)-invariant spin liquids constructed using 
fermionic spinous, the uniform spin liquid has the low- 
est variational energy in the relevant model with ring 
exchanges)^ It consists of spinous hopping on triangular 
lattice with no fluxes and thus having Fermi surface. In 
the full theory, the spinons are coupled to a U(l) gauge 
field,^i23i24^^j27 rpj^jg -g g^jjyg directly. In order 

to make progress we solve the problem in the mean field 
theory ignoring the gauge field and obtaining a system of 
free fermions hopping on the triangular lattice. We also 
go beyond mean field by using Gutzwiller projection. 

Specifically, in the mean field, we consider free fermions 
hopping on the triangular lattice in the presence of a 
missing site, a line boundary, a finite density of miss- 
ing sites, and also a random distribution of hopping am- 
plitudes. These are models of non-magnetic disorder in 
terms of what the spinons see. 

If {V'n(*)} is the set of single-particle wavefunctions, it 
is easy to show that the local susceptibility at tempera- 
ture T is given by 

X« = ^^^Y.\^n{^rf{en){l- f{en)) , (2) 

n 

where / is the Fermi function /(e) = 1/(6^*^"^'/"^ -I- 1). 
In each model of impurities, we obtain the wavefunctions 
and use this formula to obtain the local susceptibility. We 
consider various kinds of disorder in turn. We present the 
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TABLE II; Series coefficients a„ of x ~ ^ X/n=o 4"+ °(n+i)! (t)" ^'^'^ susceptibility Xuniform of tlie pure triangular lattice and 
for the five closest inequivalent neighbors Xi to a boundary as indicated in Figure |3] 



n 


Xuniform 


Xi 


X2 


X:5 


Xi 


X5 





1 


1 


1 


1 


1 


1 


1 


-12 


-8 


-12 


-12 


-12 


-12 


2 


144 


72 


120 


144 


144 


144 


3 


-1632 


-752 


-896 


-1440 


-1632 


-1632 


4 


18000 


8640 


5120 


9040 


16080 


18000 


5 


-254016 


-103488 


-108960 


-37248 


-127296 


-230976 


6 


5472096 


1497440 


3972864 


2808736 


1342432 


3429216 


7 


-109168128 


-29967872 


-58795776 


-109978368 


-46504448 


-21990912 


8 


818042112 


553745664 


-912840192 


598482432 


1331324928 


-895304448 


9 


17982044160 


-4034237440 


35460869120 


74136878080 


-6674631680 


-1900426240 


10 


778741928448 


-38283289088 


1453883081728 


-796283040256 


-765905530368 


2445141614080 


11 


-90462554542080 


-6599243882496 


-88646526167040 


-131119323998208 


1918538846208 


-58857804742656 


12 


829570427172864 


433688769173504 


-1170019148326912 


3744417183383552 


1814576120913920 


-3956791382702080 



results, a very basic discussion, and leave proper discus- 
sion of the possible connection to the k-(ET)2Cu2(CN)3 
to a later section. 

Belovif we keep the spinon hopping t fixed and vary 
the temperature, and the presented susceptibilities are in 
units of [giisY' /t- In a more systematic calculation, the 
spinon hopping amplitudes would need to be found self- 
consistently for a given spin Hamiltonian and tempera- 
ture T. In the clean system, the self-consistent t vanishes 
above some temperature of order J (e.g., in the renormal- 
ized mean field scheme this temperature is 0.75 J). When 
t becomes non-zero, this signals that the system becomes 
correlated paramagnet, and the spinon mean field is one 
attempt to capture the growing local correlations. Below 
the onset temperature, the self-consistent t quickly ap- 
proaches the zero-temperature value, and it is this regime 
that we are describing when keeping t fixed. We can esti- 
mate the spinon hopping amplitude in the renormalized 
mean field scheme as t = 3J{fj fj) w 0.5 J. The free 
fermion susceptibility on the half-filled triangular lattice 
with such t ^ lOOK would be x ~ 10 ' 10~^ emu/mol, 
which is about a factor of 2-3 larger than the experimen- 
tal values, but is reasonable given the serious approxima- 
tions in such calculations. 

We further discuss the self-consistent approach in the 
case with random bond disorder in Sec lIVP] Before that 
in the examples below, we introduce the disorder into the 
spinon problem by hand, either by simply removing the 
links to vacancy sites, or by taking randomly distributed 
bonds. 



A. Point Impurity 

In this case we find the mean field wavefunctions nu- 
merically by exact diagonalization for system sizes up to 



80 X 80. The resulting local susceptibility curves for sev- 
eral neighbors of the vacancy are in Figure |4l 
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FIG. 4: [Color online] The uniform susceptibility and the local 
susceptibilities at the first five inequivalent neighbors of the 
impurity obtained for the system of free fermions hopping on 
the triangular lattice at half-filling. The spinon hopping t is 
kept fixed and the susceptibilities are in units of {gfiB)'^/t. 

Precise curves for the susceptibility are hard to ob- 
tain at low temperatures because of the factor /(I — /) 
in Eq. which becomes increasingly sharp as T — > 0. 
Thus fewer and fewer states near the chemical potential 
contribute and eventually the results are polluted by fi- 
nite size effects. Nevertheless, we can go to sufficiently 
low temperatures with our system sizes, and the results 
in Fig. |4] for T > O.lt essentially represent the infinite- 
volume limit. We see that as we lower the temperature, 
susceptibilities at more and more neighbors become sep- 
arated from the uniform susceptibility. 

It is interesting to look at the shape of the local sus- 
ceptibility as a function of the position. At a low temper- 
ature (T = 0.054t) this is plotted in Figure [5] as obtained 
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from the exact lattice calculation. 
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FIG. 5: Local susceptibility near impurity at T = 0.074f ob- 
tained from the free fermion calculation for size 80 x 80. 



This distribution converges to a zero temperature dis- 
tribution. We can obtain some intuition about the long 
distance behavior from a calculation treating the non- 
magnetic impurity as a perturbation. The result is 



Xioc(r) -X = ^ 



cos(2fci?r + ( 



In the clean system with the boundary we can write 
all wavefunctions explicitly. The resulting susceptibilities 
are shown in Figure [B] for the first five neighbors as a 
function of temperature and in Figure [7| for fixed T = 
0.0364i as a function of the distance from the boundary. 
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FIG. 6: [Color online] The uniform susceptibility and the local 
susceptibilities at the first five inequivalent neighbors of the 
boundary obtained for the system of free fermions. x is in 
units of (gfis)^ /t. 



Here kp is the Fermi surface location where the group 
velocity points in the observation direction r, while the 
phase (j) depends on the impurity type and strength (just 
as in the case of Friedel oscillations in metals) . The calcu- 
lation leading to this result is summarized in Appendix [XI 

In the present case, the Fermi surface is roughly a circle 
with kp w [2.67, 2.72] a~^, where a is the lattice 

spacing. Taking the above expression and plotting it on 
the lattice gives a picture looking very similar to Figure[5l 
One interesting thing to notice is that there is seemingly 
much longer wavelength along the x direction than n/kp. 
This simply comes from the fact when we evaluate the 
cos(2fci?a; -|- 0) on the lattice, it picks up similar points 
at different hills of the cosine curve because the period 
n/kp ~ 1.18a is close to one lattice spacing. 

At a finite temperature, the oscillatory power law is 
cut off at the characteristic length 



ar) = nvp/{2TTT) 



(4) 



For the half-filled band on the triangular lattice, the 
Fermi velocity vp does not vary significantly with the 
direction and is « [2.86, . . . , 2.44]ia/?i. As an exam- 
ple, for T = Q.lt the correlation length is only ^ sa 4a. 

Finally, we would like to know if this distribution, with 
one impurity per 10000, can roughly give the observed 
spectral lines in the At-(ET)2Cu2(CN)3 material. The 
answer is no, and the histograms of x Si.TC still negligibly 
narrow. 




20 30 



FIG. 7: The local susceptibility as a function of the distance 
from the boundary x (line index in Fig. O at T = 0.0364t. 
Note that the vertical scale does not start at zero and that 
the X > 2 sites lie withing ±15% window around the average. 

Continuum calculation with a circular Fermi surface 
predicts at T = 



Xioc(r) - X 



X 



-M2kpr) 



cos{2kpr — 7r/4) 



y/irk 



pr 



(5) 



where Jq is the Bessel function of the first kind and r 
is the normal distance from the boundary. The asymp- 
totic form is valid also for general Fermi surface, with 
kp denoting the momentum where the group velocity is 
perpendicular to the boundary. 
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A finite temperature cuts off the power at tiie length 
scale (,{T): Roughly, the oscillatory piece is multi- 
plied by r/[^sinh(r/^)]. Indeed a six-parameter func- 
tion oo + ai cos(a2?' — a3)/[r°'* sinh(a5r)] fits the data 
like that in Fig. [7] well with 02 ~ 'ikp, 04 « —0.5, and 
0-5 ~ 1/^1 ™ agreement with our expectations. One inter- 
esting thing to notice in Fig. [7] is the apparent period of 4 
in terms of the lattice line spacing; the wavelength in the 
continuum n/kp is very accurately 4/3 of the line spac- 
ing, so the apparent period is equal to three wavelengths, 
which is an accidental commensuration effect. 

Finally, we look at the histogram of susceptibilities. 
This depends on the size of the grain - in the present 
model, the distance between boundaries. To show an ex- 
ample, we consider the system as in Fig. [H] with bound- 
aries separated by 100 lines of sites. The result is in 
Fig. [5] The grain boundary can in principle go in many 
directions or might not be straight at all. For the particu- 
lar orientation that we have chosen, the near commensu- 
ration mentioned above plays a role at low temperatures: 
For our grain size, as we lower T, the susceptibilities 
start to fail to reach the bulk value, and thus we obtain 
a double peak in the histogram (one peak from the up 
hills of the sine curve and the other from the down hills, 
cf. Fig. [7]). This starts to happen for temperatures just 
slightly below the smallest one shown in Fig. [8l 

We do not know if this type of disorder is realistic 
in the k-(ET)2Cu2(CN)3. The chosen separation of 100 
spacings between boundaries is a rather significant disor- 
der: 2% of the sites are right next to the boundaries and 
several times more are in the immediate vicinity. Still, 
at temperature above O.K the histograms are very nar- 
row is still smaller than 5 lattice spacings), which is 
why they are shown only below this temperature. Even 
at the lowest temperature the linewidth is small, despite 
the slow decay of Xioc('') away from the boundary. Think- 
ing about the k-(ET)2Cu2(CN)3, it appears that we need 
more disorder than this and more spread across the sys- 
tem. 
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FIG. 8: Histogram of the susceptibilities for the system with 
grain boundaries separated by 100 lines of sites. 



C. 5% of Impurities 

In this section we calculate the susceptibility histogram 
for the samples with 5% of spin vacancies. It seems 
unlikely that this type of disorder is present in the k- 
(ET)2Cu2(CN)3in the form of missing ET dimers. How- 
ever, this could be a crude spin model if the electron 
charge distribution is inhomogeneous. There are other 
likely sources of disorder and this case represents a situ- 
ation when the disorder is point-like. The resulting his- 
tograms are in Figure[9l We see that peak is quite broad, 
more like the experimental curves, and it spreads as we 
lower the temperature. 
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FIG. 9: Histogram of the susceptibilities for the system with 
5% of vacancies. Note especially that the peaks broaden as 
we lower the temperature. The curves are obtained for 60 x 60 
samples and are in the thermodynamic limit for temperatures 
shown. 



D. Random Spinon Hopping Amplitudes 

In this section we take a model of non-magnetic dis- 
order where the spinon hopping amplitudes are random 
and uniformly distributed in the range [t — A, i 4- A]. 
The resulting histograms for A/t = 0.2 are in Figure [TUl 
Note specifically that the peak of the histogram does not 
change much below roughly T ~ 0.44t (the peak does not 
spread). 

We also tried to make disorder more point-like and 
see if this would cause the peak to spread, as it did for 
the point-like missing sites above. We find that this is 
indeed the case for the following simple choice: Take 95% 
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of bonds to have one value and 5% to have twice as large 
value (results not shown). 

The fact that point-like disorder spreads the peak upon 
lowering temperature can be understood as follows. As 
we lower T the correlation length ^(T) grows and more 
and more sites start to "feel" the impurities and have 
susceptibility substantially different from the bulk value. 
This will be happening until the correlation length be- 
comes somewhat larger than the typical spacing between 
the point impurities. For the samples studied we can 
go as low as T = 0.05i (before finite size effects set in), 
and at this T the correlation length is roughly ten lattice 
spacings. Thus we don't expect the histogram in Fig. [9] 
with 5% vacancies to broaden much as the temperature 
is lowered further. 
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FIG. 10; Histogram of the susceptibilities for the system with 
random values of bonds with randomness ±20%. Note es- 
pecially that the peaks do not broaden significantly below 
T ~ OAit. The curves are obtained for 60 x 60 samples and 
are in the thermodynamic limit for temperatures shown. 

So far we took a specific fixed distribution of spinon 
hopping amplitudes and calculated susceptibilities. As 
we argue below, this is reasonable for a model of non- 
magnetic disorder where couplings in the spin Hamilto- 
nian have some randomness (assuming it is not too strong 
and the spin liquid state remains stable). To treat the 
system more properly, the spinon hopping amplitudes 
should be calculated self consistently. For example, in 
the Heisenberg model with exchanges Jrr' , popular mean 
field self-consistency conditions read i*^, = Jrr'ifrfr'), 
the bond expectation value is for one spin species and 
Jrr' is proportional to Jrr'-, e.g., in the so-called renor- 
malized mean field scheme one takes J = 3J. In the 
clean system, the above self-consistency condition has 



non-trivial solution once the temperature is lower than 
T = J/4. Rather quickly below this temperature the 
spinon hopping becomes large and comparable to the 
zero-temperature limit. However, in this specific treat- 
ment, the spin liquid with the Fermi surface is not a sta- 
ble solution and other spin liquids perform better. Fur- 
thermore, as is known from early slave particle studies, 
the best states in this mean field arc dimcrized. In par- 
ticular, if we try to solve the self-consistency equations 
by iteration starting from a random initial trr' , these run 
away towards some dimerized solutions. 

As discussed in the introduction, we expect that there 
are additional spin interactions such as ring exchanges 
that stabilize the zero flux spin liquid against other spin 
liquid states and dimerized states. Ref. H presented a 
schematic mean field argument how the ring exchanges 
achieve this. The self-consistency condition is modified 
to 



= Jrr' iflfr') + J2 ^^^P^frfs) {flfs'){flfr') , (6) 



where s, s' are all the sites so that each expectation value 
is for two sites on a bond, and such sum effectively covers 
all four-site rhombi P ~ [rss'r'] on the triangular lattice 
that contain the bond rr' . The couplings Kp are propor- 
tional to the ring exchanges acting around the rhombi. 

Ref. applied the above scheme to the clean system 
at T = and found that the uniform spin liquid with no 
fiuxcs in the hoppings is a stable solution for K / J > 9.9. 
Note that the parameters J and K arc related to the 
microscopic Heisenberg and ring exchanges by disparate 
numerical factors, and Ref. @ contains more details in 
what sense such I'C / J values are reasonable in the study 
of the spin liquid. Here we mainly use this scheme to have 
a starting point where the uniform state is stable in the 
clean system and see crudely how the randomness in the 
microscopic parameters like J translates to randomness 
in the spinon hopping amplitudes. 

To get some understanding of the self-consistent distri- 
bution of trr' and its temperature dependence we iterated 
Eq. ^ until convergence for the following system pa- 
rameters. We took Kp = 15 everywhere and took a uni- 
formly distributed Jrr' from the interval [1 — A, 1-1- A] for 
two values A = 0.05, 0.2. We find that once the nontrivial 
solutions appear, which happens quickly below T ^ 1/4 
(in units of J), the distribution of trr' is essentially in- 
dependent of the temperature and has the same width 
in relative terms as the distribution of Jrr' but is a bit 
more rounded. Calculating histogram of susceptibilities 
for this distribution gives roughly the same result as cal- 
culating it for the box distribution of trr' 's of the same 
relative width as that of Jrr'- This provides some jus- 
tification to the preceding models of disorder where we 
simply put randomness into the spinon problem by hand. 
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E. Gutzwiller Wavefunction Study of the Local 
Magnetization 

Let us discuss the spin liquid picture beyond the 
mean field. One way to proceed is to consider ef- 
fective gauge theory description where spinons interact 
with the emergent gauge field. One expects that the 
power laws in Eqs. pisp arc modified by the gauge field 
fluctuationSf^i^i^ but reliable quantitative information 
is lacking. 

In this work we go beyond the mean field by Gutzwiller 
projection. In the ID case, this essentially reproduces 
exact result^ for the Hciscnbcrg system with a nonmag- 
netic impurity. However, in 2D the Gutzwiller projection 
alone likely does not capture all important fluctuations in 
the low-energy theory^ Nevertheless, by working directly 
with the physical spin variables, it gives quantitatively 
more plausible results than the mean field. 

Specifically, we consider local magnetization distribu- 
tion in a partially polarized state both in the mean field 
and after the projection. We used this approach to 
study non-magnetic impurities in a kagome spin liquid 
in Ref. [s^ (this reference also contains more discussion 
on the connection to the local susceptibilities). 

It is well known that the Gutzwillcr-projccted Fermi 
sea is an excellent trial wavefunction for the ID spin- 
1/2 chain, and we test our approach in this case. Fig- 
ure [TT] shows results for a chain with open boundaries. 
In the mean field, the local magnetization is xioc(a;) ~ 
1 — cos{2kpx) = 1 — (—1)^. The projection dramatically 
enhances the staggered component. In the Heisenberg 
chain, Eggert and Affleck^^ predict that the staggered 
component in ^loc grows as ^/x away from the boundary 
at T = 0, and the Gutzwiller-projected state appears to 
capture this result in the toioc- This dramatic behavior 
of the xioc near a non-magnetic impurity has been used 
to explain broad lines in spin- 1/2 chain compounds even 
with small density of impurities. "^^ 

Our initial hope was that the 2D spin liquid, which is 
also a projected Fermi sea state and whose full theory 
shows enhanced spin correlations at 2fci?, could similarly 
produce broad xioc histograms around small density of 
impurities. However, it appears that quantitative aspects 
in 2D are such that small impurity concentration docs not 
give large line broadening. 

Specifically, in the 2D spin liquid case, we studied both 
a single impurity and a line boundary and found that 
™ioc — m is enhanced by projection only by a fixed nu- 
merical factor of about two. Figure [T^ shows representa- 
tive results in the line boundary case, where the impurity 
perturbation is the largest of the two cases, cf. Eqs. ([3]) 
and ([5]) . The triangular lattice is constructed by stacking 
24 chains of length 24, with periodic boundary conditions 
along the chains and open boundaries in the perpendicu- 
lar direction. In an effort to bring out more effect, the ex- 
cess spin-up population occupies orbitals near the patch 
where the Fermi velocity is normal to the boundaries, 
since it is this fc-spacc region that is responsible for the 
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FIG. 11: [Color Online] Local magnetization profile in a ID 
chain of length L — 129 with open boundaries before and af- 
ter Gutzwiller projection of a partially polarized spinon Fermi 
sea state, here with A*'f = 65, — 64. The excess up spin 
occupies orbital sin(fcF2;) = sin(7ra;/2), producing a staggered 
component in xioc with a constant amplitude. The projec- 
tion dramatically enhances the staggered component, which 
now grows as ^/x away from the boundary, in agreement with 
Ref. for the Heisenberg chain. Such effects can produce 
broad xioc histogram even for small impurity density in the 
ID chain. We want to contrast this with the 2D case, where 
we find only a fixed numerical enhancement. 



power law in Eq. ([5|). Note that because of this spe- 
cial population, the mean field amplitude of oscillations 
is larger compared to the case with thermal population 
of orbitals (irrespective of orientation) in Fig. [T] Never- 
theless, the figure shows that the Gutzwiller projection 
gives only a fixed enhancement over the mean field by a 
factor of about two. Different system sizes and orbital 
populations do not change this result qualitatively. Sim- 
ilar numerical enhancement was observed in our Kagome 
study^ where we also argued for it using renormalized 
mean field thinking. 

To conclude, we expect that all our mean field results 
will experience a similar numerical enhancement in the 
Xioc — X by the projection, so the histograms will be 
broader by about a factor of two. In particular, we can 
find negative local susceptibilities, despite the mean field 
giving only non-negative xioc- But unlike the ID, we are 
not able to get small density of impurities to produce 
broad lincshapes. One caveat here is that, as wc have 
already mentioned, the Gutzwiller projection in 2D does 
not capture the full gauge theory, and it could be that 
the effects of gauge fluctuations are much more dramatic 
(e.g., see footnote^). This could happen if the actual 
spin liquid phase has much stronger correlations than 
the mean field prediction, but at the moment we do not 
know how to address this better quantitatively. 
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FIG. 12: [Color Online] Local magnetization profile in a trian- 
gular lattice with grain boundaries before and after Gutzwiller 
projection of a partially polarized spinon Fermi sea state. The 
lattice is constructed of 24 chains of length 24, with periodic 
boundary conditions along the chains and open BC in the 
perpendicular direction, see the drawing in Fig. [U] The par- 
tial spin polarization is obtained starting from an unpolarized 
state and depopulating 5 spin-down orbitals and occupying 5 
spin-up orbitals near the Fermi patch whose normal is per- 
pendicular to the boundary. 



V. INHOMOGENEOUS KNIGHT SHIFTS IN 
THE METALLIC PHASE 

In Sec. llVi we used free fermions as a mean field for the 
spinons. Assuming the spin liquid is appropriate in the 
insulator, the treatment further neglects gauge fluctua- 
tions and is a crude approximation that can change qual- 
itative long-distance behavior. However, this was best 
we could do to get some quantitative estimates of local 
properties. 

In fact, the free fermion analysis applies more read- 
ily to the metaUic phase of the k-(ET)2Cu2(CN)3. The 
fermions are now electrons themselves, and mean field 
is reasonable in the Fermi liquid regime. We use the 
same formula Eq. ([2]) to calculate the local susceptibil- 
ities. Also, the analytical results Eq. ^ for the long- 
distance behavior of Xioc('') away from a single impurity 
and Eq. ([5|) away from a boundary hold in the metal. 

Here we are interested in modelling weak disorder in 
the metallic phase and connecting with the ^^C NMR 
measurements. In this case, we calculate the histogram 
of susceptibilities in the presence of random on-site po- 
tentials taken to be uniformly distributed in an interval 
[-W, W]. The result for W = 0.3ie is in Figure [H 

Here te is the electron hopping amplitude and is in 
excess of 50 meV. This model of disorder is reasonable 
in the following sense. Note that the disorder strength 
W should be compared to the band width, which is sev- 
eral times te, so W = 0.3te is a rather weak disorder. 
A crude Born approximation of the electron mean free 
path due to elastic scattering in the 2D system gives 
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FIG. 13: Histogram of the susceptibilities in the metallic 
phase, where we model disorder by random on-site potentials 
uniformly distributed in an interval [—W, W\ with W = O.Ste- 
The curves are obtained for 60 x 60 samples and are in the 
thermodynamic limit for temperatures shown. 



kpl w h^Vp/ {A6fi'^), where A is the area per site and Sfj!^ 
is the variance of the on-site potential. For the triangu- 
lar lattice at half-filling we have kpl « St'^/S^'^, which 
gives kpl ~ 250 for the above disorder. Direct numerical 
estimate of the lifetime of momentum eigenstates in the 
lattice system gives comparable number. Residual resis- 
tivities in the k-(ET)2Cu2(CN)3 in the metallic phase 
imply kpl > 50 and larger, so the disorder that we use is 
reasonable. 

Examining the local susceptibility histograms in 
Fig. [T31 we see that below roughly T = OMe the peak 
no longer spreads. This T is comparable to the room 
temperature, and indeed the available data in the metal- 
lic phase shows little temperature dependence of the 
linewidth»^ Our model linewidths are reasonable, even 
though we cannot read off reliably the inhomogeneous 
broadening component from the lines in Ref. |^. 

In the free fermion case of spinons in the previous sec- 
tion we found that making the bond disorder more point- 
like, by changing value of only a fraction of bonds from 
the uniform value, the peak of the histogram spreads as 
we lower the temperature over a wider temperature range 
compared to the case where some randomness is present 
on all bonds. In this section we also studied whether sim- 
ilar effect takes place for electrons with random on-site 
potentials. We set 95% of chemical potentials to zero and 
5% to a positive value, chosen to be OM in one run and 
2t in the other. We found that indeed the peaks spread in 
this case too. This reinforces the conclusion made above 
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that more point-like disorder causes the xioc histograms 
to have stronger temperature dependence of the spread. 



VI. DISCUSSION 

We summarize the main results with an eye to connect 
with the NMR experiments in the k-(ET)2Cu2(CN)3. 

First, our high-tcmpcrature series study shows that the 
local susceptibility near non-magnetic impurities such as 
vacancies or grain boundaries can deviate sizably from 
the bulk value. We studied specifically the Heisenberg 
model with vacancies and obtained quantitatively accu- 
rate results down to T ^ J/4. Even at this low tem- 
perature, the local susceptibility is modified perceptibly 
only within few lattice spacings of the defects, cf. Figs. [T] 
and m On the other hand, the ^"^C NMR experiments 
show broadening that develops gradually from high tem- 
peratures; the linewidth roughly doubles going from T ^ 
J ~ 250 K down to 50 K, at which the FWHM aheady 
corresponds to about 20-30% of the bulk susceptibility. 
In our study, the few close neighbors of the defects show 
similar deviation in xioc- However, for a small density 
of defects, which was our initial assumption, essentially 
all sites would be sufficiently away from impurities and 
we would not be able to obtain comparably broad his- 
tograms. 

We studied vacancies as a model of local nonmag- 
netic disorder, but we do not expect significant changes 
for other types of disorder such as random bonds. The 
Heisenberg model is also not entirely adequate at T = 0, 
where multi-spin exchanges likely affect the ground state, 
but the role of such terms in a nominally spin model is 
less important at higher temperatures. The Heisenberg 
model is thus a reasonable choice at such temperatures 
and was already used successfully to understand the bulk 
spin susceptibilityFiiii From the high-temperature study 
with defects, we are led to ask if there is more disorder in 
the system than originally thought. To be able to repro- 
duce the T = 50 K lines, we seem to need the disorder 
strength comparable to that of few percent vacancy con- 
centration. Unfortunately, the way the high-temperature 
series work, we cannot study directly more realistic mod- 
els of disorder such as random bonds, but our work with 
vacancies gives a rough idea. 

Next, we considered the spin liquid with spinon Fermi 
surface as a plausible description of the correlated param- 
agnet in the temperature range below 50-100 K and down 
to few Kelvin. This is a serious assumption, and even 
within it we can do quantitative calculations only in the 
mean field approximation, supplemented by Gutzwillcr 
renormalizations. Proceeding nevertheless, in such spin 
liquid at low temperatures, xioc('') — X decays with slow 
power laws away from defects, cf. mean field Eqs. ([3])^ 
and (O, and many sites can be potentially affected by 
impurities. However, quantitative aspects appear to be 
such that we would not get visibly broadened histograms 
even at T = unless there is a significant density of 



impurities. If we postulate a sizable disorder, we can 
get Xioc histograms comparable to the experimental ones 
in the temperature range 50 to 10 K. We find that the 
variation with temperature depends on the type of dis- 
order. If the disorder is uniformly spread, e.g., all bonds 
are random in some range, the peak stops broadening 
at a relatively high temperature of about half the over- 
all spinon hopping amplitude. On the other hand, for a 
more point-like disorder, the peak keeps broadening to 
much lower temperatures. Thus, if the disorder strength 
is fixed, in order to get significant temperature depen- 
dence of the linewidth in the spinon analysis we would 
need a point-like disorder. 

The free fermion mean field applies directly to 
the metallic side of the phase diagram of the k- 
(ET)2Cu2(CN)3, which appears for pressures above 
0.4 Gpa. The mean field fermions are now the elec- 
trons themselves. What is observed^ are essentially 
temperature-independent ^"^C NMR lines. As a reason- 
able type of disorder in this case we took a random distri- 
bution of the chemical potentials. Specifically, for a box 
distribution [— O.S^e, O.S^e] , which gives a reasonably large 
kpl ^ 250, we find sensible and temperature-independent 
Xioc histograms. On the other hand, making the disorder 
more point-like, wc find that the histograms broaden as 
we lower the temperature. Thus on the metallic side the 
NMR lines suggest a uniformly spread disorder. 

Returning to the Mott insulator side, it appears that 
there is more disorder here than in the same system on 
the metallic side. Furthermore, if the disorder is fixed 
as is reasonable in the metal, the metallic side suggests 
it is uniform and not producing the broadening of the 
lines, and so this should also be the case on the Mott 
insulator side, which contradicts the experiments. How- 
ever, let us assume for a moment that the disorder is 
point-like. The peak in Fig. [9] produced for 5% vacancies 
broadens by about a factor of two or three in the range 
T = OM to Omt. Using t - 100 K, in the experiment 
this correspond to the range 50 K to 5 K where we see 
broadening by about a factor of three which is thus a rea- 
sonable agreement. As mentioned, we don't expect much 
broadening below O.OSt because the correlation length is 
already about ten lattice spacings which is comparable to 
the distance between impurities in this example. How- 
ever the experiment shows very strong broadening be- 
yond that. It might be that the origin of the broadening 
is not disorder and that our spin liquid picture is not 
adequate and some other phase is emerging at low tem- 
peratures, perhaps with incipient magnetic order F^^i^i^^ 
There is a way to explain the observed phenomena within 
spin liquid picture, but it seems to require that the disor- 
der strength effectively grows as the temperature is low- 
ered. Below we speculate on how this may come about, 
but more experimental input is needed. 

One source of disorder mentioned in the literature 
for the (ET)-based organic superconductors is ethy- 
lene group disorder J^iii This was particularly discussed 
for the K-(ET)2Cu[N(CN)2]Br material, where signifi- 
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cant sample to sample variations and cooling rate de- 
pendence were observed. However, recent studie a^^i"'^^ 
suggest that the amount of such disorder in the k- 
(ET)2Cu[N(CN)2]Br is small at low temperatures and 
that perhaps the insulating polymeric layer is involvedf^ 

For the spin hquid material k-(ET)2Cu2(CN)3, 
literaturo^i^i^ mentions that one of the (CN)~ groups in 
each unit cell is orientationally disordered and that such 
structural disorder can generate random electrostatic po- 
tential throughout the latticci^ If such disorder is in- 
deed involved, the question is then why it does not have 
comparable pronounced effects on the metallic side. One 
possible explanation is good screening of charged impu- 
rities in the metal and progressively weaker and eventu- 
ally absent screening in the Mott insulator. For exam- 
ple, in the metallic phase, the Thomas-Fermi screening 
length is small, Xtf = I / ^/ATiehAep) 1 A, where 
we estimated the density of states at the Fermi level by 
v{eF) = 0.28/(tew), te = 50 meV is the electron hopping 
amplitude, and v = 850.6 A'^ is the 3D volume per tri- 
angular lattice site. On the insulator side, an accurate 
calculation is harder to make. Using semiconductor lan- 
guage, we can estimate A = ^ ksT / (Aire^n) , where n is 
the number of thermally excited charge carriers. The re- 
sistivity of the insulator at ambient pressure increases by 
about four orders of magnitude when the temperature is 
decreased from room down to 25 K. Taking this as an or- 
der of magnitude measure of the change in the density of 
carriers, we get A ~ 10—100 A, which is about several lat- 
tice spacings. So one scenario is that the charged disorder 
is still well screened at room temperature but gradually 
becomes more visible below 100 to 50 K, with the screen- 
ing essentially absent below about 10 to 5 K. It might be 
that this type of disorder is more point like, which further 
enhances the broadening of the xioc histogram on the in- 
sulator side, but is screened on the metallic side where 
only weak and more uniformly spread disorder remains 
that does not cause significant histogram broadening at 
lower temperatures. In Appendix [B] we briefly discuss 
how the charged disorder in the electronic system may 
translate to that in the spin model for the insulator. 

At present, we do not know how to estimate the 
strength of disorder in the system and whether the above 
scenario is reasonable. Unlike the metallic phase, we can- 
not use the electrical resistivity as a measure of impu- 
rity scattering. We want to remark though that if the 
disorder is not too strong so that the spin- 1/2 model 
with say random couplings is applicable, our spin liq- 
uid construction is still a viable candidate for the ground 
state. As discussed in Sec. IIVDI we can accommodate 
the bond disorder by adjusting spinon hopping ampli- 
tudes, which are now non- uniform. From the point of 
view of spinous, they arc now scattered by this disorder. 
Interestingly, if the corresponding elastic mean free path 
is sufficiently small, the thermodynamics of the spinon- 
gauge system with such diffusive spinous differs from the 
clean casc i^^i^^ For example, the specific heat behaves as 
C'diffusivo ~ Tlog{l/T) as opposed to the clean system 



C'cican ~ T^/'^, while the thermal conductivity behaves as 
^diffusive ^ r^/^ as opposed to Kcioan ^ T^/'^. The ther- 
mal conductivity measurements could potentially give an 
independent estimate of the degree of disorder in the in- 
sulator, while at present the ^'^C NMR is our only window 
onto the disorder. 
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APPENDIX A: PERTURBATIVE CALCULATION 
IN IMPURITY STRENGTH 

For a weak non-magnetic perturbation we find, to lin- 
ear order. 



Xioc('7) 



(Al) 



^ /(efc)[l~/(efc)]-/(6fc+,)[l-/(efc+,)] ^^^^ 
where e(fc) describes the clean system band structure, 



while 



''imp 



is the appropriate impurity matrix element. 



for which we neglect any dependence on the momenta. 
Analyzing the integral in the T — > limit, the q- 
depcndence has singularities on the 2k p surface of the 
form 



e{q > 2kF) 



2TTvl^Jc{q-2kF) 



(A3) 



Here is the Fermi velocity and c the curvature at 
the appropriate Fermi surface patch. Going back to real 
space, we obtain Eq. ^ with A = — Ui,„p/(27r^i;|,c) and 
= 0. 

The 2kF oscillations of the local susceptibility of free 
fcrmions in the presence of impurities can also be related 
non-perturbatively to the calculation of the Friedel oscil- 
lations in the local density ; '^^i'^^ 



n 

dp\oc{i) dp\oc{i) 



Xioc(i) 



djjL 



dkF 



(A4) 
(A5) 



This gives one power of r slower decay at large distances 
in xioc('') than in p\oc{i') away from defects. 
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APPENDIX B: FROM RANDOM ELECTRON 
POTENTIAL TO RANDOM HEISENBERG J 

Here we discuss schematically how disorder in the mi- 
croscopic electronic model translates to disorder in the 
effective spin description of the insulator. For an illus- 
tration, consider a Hubbard model 

H = U '^nrijTlri-^trr'cl^Cr'cr+'^Vrcl^Cra- , (Bl) 
r rr' r 

with hopping amplitudes trr' and on-site potentials iv- 
Large U opens a charge gap. To leading order in 1/U, 
we obtain spin-1/2 model with Heisenberg couplings 

4t2 , 

= tjU )1 ■ ^^^^ 

Thus, if we have either random hopping amplitudes or 
random potentials, the Heisenberg couplings are also ran- 
dom. 



The randomness in the electron hopping amplitudes 
translates directly to randomness in the spin exchanges 
- the latter is even larger in relative terms. On the other 
hand, if the disorder is in the potentials and if these 
are much smaller than [/, this randomness is effectively 
renormalized down. This is natural since the variables in 
the spin model are charge-neutral and should be oblivious 
to such potential disorder. Still, the random potentials 
are seen via virtual excitations and induce non-magnetic 
randomness in the spin system. 

If we need to include multi-spin exchanges so as to sta- 
bilize the spin liquid, there will be some randomness in 
these as well. Despite the randomness in the effective spin 
model and irrespective of its electronic origin, as long as 
it is not too strong, we can proceed with the spin liq- 
uid construction but now the spinon hopping amplitudes 
become non-uniform. This was discussed in Sec. IIVDI 

Finally, if the random potentials become comparable 
to the Hubbard U and lead to significant variation in 
the electron density, the very spin model thinking can 
become inappropriate. 
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